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Electroactive soft elastomers require huge electric field for a meaningful actuation. We demon- 
strate that this can be dramatically reduced and giant deformations can be produced by application 
of suitably chosen heterogeneous actuators. The mechanism by which the enhancement is achieved 
is described and illustrated with the aid of both idealized and periodic models. 

PACS numbers: May be entered using the \pacs{#l} command. 



Electroactive polymers (EAP) are capable of large de- 
formations in response to electric stimulus. A sketch of 
a planar actuator is shown in Fig. [l] The top and bot- 
tom faces of the soft dielectric are covered with compli- 
ant electrodes [1 inducing an electric field through the 
material. The resulting Maxwell stress leads to the de- 
formation of the material. The variety of possible ap- 
plications of these "artificial muscles" motivated an in- 
tensive search for appropriate polymers. Indeed, recent 
experimental studies achieved remarkable milestones in 
terms of the magnitudes of the actuation strains [IHZ]- 
In parallel, the concept of enhancing the responsiveness 
of EAP devices by means of snap-through unstable mech- 
anisms was examined too [SHIO]. However, these studies 
did not tackle the main limitation of EAPs, namely the 
huge electric fields needed for meaningful actuations. 

We address this challenge and investigate a mechanism 
by which the exciting electric field can be reduced by an 
order of magnitude. In this regard we recall that re- 
cent experimental works f7\ [Tl involving soft elastomers 
with high dielectric particles demonstrate an improved 
response of the heterogeneous systems. In agreement, 
theoretical studies [12^ 13^ of idealized heterogeneous di- 
electrics predicted an enhancement of the electromechan- 
ical coupling. In this work the non-linear theory of elec- 
troelasticity at finite strains p!4HT7] is adopted and an 
exact analytical solution for an idealized heterogeneous 
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FIG. 1. A sketch of a planar EAP actuator. 
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FIG. 2. A sketch of the idealized layered model. 



system is deduced. This, in turn, sheds light on the mech- 
anism that leads to the improved coupling and motivates 
an investigation of more realistic microstructures. By ap- 
plication of the finite element (EE) method correspond- 
ing periodic models are examined and the improvement 
in the electromechanical coupling is quantified. 

The deformation of the material is characterized by 
the deformation gradient Fij = dxi/dXj^ where Xi and 
Xj are the position vectors of a material point in the 
deformed and reference states, respectively. The electric 
field at a point is Ei = —d(j)/dxi^ where <p is the electric 
potential. In incompressible and isotropic neo-Hookean 
dielectrics the electric displacement and the total stress 
tensor are 



eocEi and crij = jiFikFjk ^ eQeEiEj - p5ij, (1) 



where eo is the vacuum permeability, e is the dielectric 
constant, ji is the shear modulus, p is the pressure, and 5ij 
is the Kronecker delta. In the expression for the stress the 
first term is the mechanical stress and the electrostatic 
Maxwell stress tensor is the second term. Assuming a 
quasistatic deformation, no magnetic fields and no body 
forces, the governing equations are 



—— = and —-^ 

OXj ox^ 



0. 



(2) 



In the absence of free charges at the interfaces the electric 
field continuity conditions are \Ei\ rhi = and {Dij hi = 
0, where [•] = (•)+ — (•)" is the jump across the in- 
terface and hi and rhi are the unit vectors normal and 
tangent to the interface, respectively. The corresponding 



mechanical continuity conditions are {Fij 



and 



an rii 



0. 
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FIG. 3. A sketch of the periodic model. 



Figure [2] depicts an idealized layered material with al- 
ternating soft isotropic and anisotropic layers (phases 2 
and 1). The anisotropic layers are themselves laminated 
structures made out of alternating "sublayers" of a stiff 
material with high dielectric modulus and a soft mate- 
rial with low dielectric modulus. Four parameters char- 
acterize this microstructure: the lamination angle with 
respect to the electrodes plane and the volume fraction 
of the isotropic layers (62 and 0^2), and the lamina- 
tion angle and volume fraction of the stiff sublayers in 
the anisotropic layer (0i and ai). It is assumed that 
the thickness of the sublayers is an order of magnitude 
smaller than the thickness of the enclosing layer. This 
"scale separation" allows to solve analytically the associ- 
ated coupled homogenization problem [HI [19], and to de- 
termine the planar mechanical response to electrostatic 
excitation Eq = (/>o/<^, where is the electric poten- 
tial between the electrodes and d is the distance between 
them in the undeformed state. The parameters of the two 
models examined in this work are summarized in Table [H 
Model II follows an optimization process in the limit of 
infinitesimal deformation that resulted in an extremely 
thin isotropic layers [13]. Model 12, with thicker isotropic 
layers, is reminiscent of the periodic models that were in- 
spired by the predictions of the optimized one. 

A finite element model of a particulate periodic mi- 
crostructure is shown in Fig. [3j The model consists of 
thin isotropic layers at an angle O2 between bulky stacks 
of elongated stiffer inclusions that are tilted at an angle 
9i. The highlighted inclined rectangle is a representa- 
tive unit cell (RUG) of the model P5 with 5 inclusions. 
A denser model, PIO, with 10 inclusions in the RUG is 
analyzed too. The aspect ratios a/b are 16.35 and 32.7 
for the models P5 and PIO, respectively. The simula- 
tions were accomplished with appropriate in-plane pe- 
riodic displacement and potential boundary conditions. 



TABLE I. Microstructure parameters 

Model ai 0L2 Oi O2 

"Tl 0^416 0008 63l^ 27^5^ 

12 0.6 0.06 63.1° 27.5° 
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FIG. 4. Principal stretch ratio versus the electrostatic ex- 
citation. The black and magenta curves correspond to the 
idealized models II and 12. The red and blue curves are the 
results of the simulations with models P5 and PIO. The green 
curve is the response of a homogeneous VHB-4910 actuator. 

The non-linear coupled problem was solved by applica- 
tion of GOMSOL FE code. 

Properties of the widely used 3M VHB-4910 scotch 
[6 were chosen for the soft phase in both, the idealized 
and the periodic, models. The properties for the stiffer 
high dielectric inclusions are comparable with those of 
Polyaniline (PANI) [11 . The values of the physical con- 
stants are listed in Table (Til 

Predictions for the principal actuation stretch as func- 
tions of the applied electric excitation are shown in Fig. [4] 
The response of the homogeneous VHB-4910 actuator is 
determined via the relation [19 

X={l-eeoEi/^i)-'^\ (3) 

The optimized microstructure II displays the best re- 
sponse. Nonetheless, even the more realistic periodic mi- 
crostructures exhibit a giant enhancement of the response 
relative to that of the homogeneous actuator. Specifi- 
cally, at 17MV/m the periodic model PIO attains a 20% 
strain in comparison with the 2% strain of the homoge- 
neous actuator. 

The amplification of the electromechanical coupling is 
a tricky two-parts mechanism that highlights the roles 
of the fluctuations in the electric field and anisotropy. 
Firstly, due to the continuity of the normal component 
of the electric displacement, the ratio between the mean 
electric fields in the anisotropic and the isotropic layers is 
inverse proportional to the ratio between their dielectric 
moduli. Gonsequently, the electric field in the anisotropic 



TABLE II. Material constants 

Phase e /x [MPa] 

VHB 4910 6^5 0^2 

Polyaniline (PANI) 6500 2700 
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FIG. 5. The mean electric fields in the isotropic layers versus 
the electrostatic excitation. The black and magenta curves 
correspond to the idealized models II and 12. The blue and 
red curves are the results of the FE models PIO and P5, and 
the error bars represent the standard deviation in the field. 
The green curve depicts the field in a homogeneous VHB-4910 
actuator. 

layer is small, and since the volume average of the electric 
field is equal to Eq^ the magnitude of the electric field in 
the isotropic layers is proportional to the applied field 
divided by their volume fraction [13 . In Fig. [s] the mean 
electric fields in the isotropic layers are shown. Indeed, 



the intensities of these fields are approximately Eo/a^'^\ 
As the number of stiff inclusions in the RUG increases, 
the curves for the periodic models approach the one for 
the idealized model 12. 

The large electric field results in a large electrostatic 
stress that tends to stretch the isotropic layer in the 
transverse direction in a manner reminiscent of the one 
shown in Fig. [T] Roughly speaking, due to the large fiuc- 
t nations in the electric field, the isotropic layers act like 
"micro-actuators" . This brings us to the second stage of 
the amplification mechanism. Since the anisotropic lay- 
ers are made out of alternating stiff and compliant sub- 
layers, their compliant mode corresponds to a shear of 
the soft sublayers with rotation of the stiffer ones. This 
mode amounts to an extension at 45° to the sublayers 
plane, which is quite close to the angle O^^^ — O^^^. Thus, 
when the micro-actuators deform due to the intensive lo- 
cal electric field, they stretch the bulky anisotropic layers 
along their soft mode, and the actuator expands in a di- 
rection transverse to that of the electric excitation. 

The amplification mechanism is illustrated in Fig. [6| 
where the intensification of the electric field in the micro- 
actuators, their resulting deformation perpendicular to 
the direction of the electric excitation, and the rotation 
of the elongated stiff inclusions can be appreciated. We 
conclude noting that the proposed mechanism exhibits a 
ten-fold enhancement of the electromechanical coupling, 
thus providing a method to overcome the primary obsta- 
cle in the development of EAP actuators. 
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